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We show that saturation can manifest itself in central-forward dijet 
production in p-A collisions. In spite of large transverse momenta of the 
jets, the almost back-to-back dijet configurations are able to probe gluon 
density at low x and low kt ■ We perform our study in the framework of high 
energy factorization with the unintegrated gluon density given by a nonlin- 
ear QCD evolution equation. We show that the formalism can successfully 
account for features measured in e-p and p-p data and it predicts significant 
suppression of the central-forward jet decorrelations in p-Pb compared to 
p-p, which we attribute to saturation of gluon density in the nucleus. 

PACS numbers: 13.85.Hd 



1. Introduction 

Nonlinear effects of QCD dynamics are expected to become increasingly 
important as one goes for a more and more forward production. This is 
because the forward region corresponds to the density of the incoming gluon 
being probed at small values of the longitudinal momentum fraction, x, for 
which a typical kt of the gluon is of the order of, or smaller than, the 
saturation scale Qs{x). To further enhance the relative importance of the 
saturation region, one may go to the proton- nucleus (p-A) collision as the 
saturation scale in the nucleus is expected to be significantly higher than in 
the proton. 

The LHC experiments were designed predominantly to measure hard 
final states like high-pt jets or vector bosons. The framework that allows 
one to study dense gluon systems by the means of looking at hard objects is 
provided by the high energy factorization formalism [T] . In this contribution 
we show that saturation effects may indeed be observed by looking at the 
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Fig. 1. Jet production in the forward region in hadron-hadron collision. 

azimuthal correlations between two leading jets from p-A collisions, one of 
which is produced in the forward direction. Such a final state probes the 
parton density of the nucleus at low x, while that of the proton at relatively 
large longitudinal momentum fraction. 



2. High energy factorization 

The leading order contribution to dijet production comes from the 2—7-2 
partonic process 

a{ki) + b{k2) ^ c{pi) + d{p2) . (1) 

The fractions of the longitudinal momenta of the initial state partons are 
related to the transverse momenta and rapidities of the final state partons 
by xi^2 = {ptiG^^^ + Pt2e^'^^) with S being the squared energy of the 
incoming hadrons. In the case of forward-central or forward-forward jet 
production, these fractions of momenta are highly asymmetric, xi ~ 1 and 
X2 ^ 1, we can therefore neglect the transverse momentum of the parton a. 
Using this fact together with the Sudakov decomposition of 4-momenta of 
the incoming partons, we arrive at the following formula for dijet production 
cross section in the high energy factorization approach 

da 

dyidy2dpitdp2tdA(j) 

a,c,d 

where k"^ = p^^ + Pf2 + '^PtiPt2 cos A(/» and A(/> is the azimuthal distance 
between the two outgoing partons. All the notation in Eq. ([2]) corresponds 
to that in Fig. [TJ Magged matrix element for the 2—7-2 process with 
one off-shell initial state gluon was calculated in [2]. Since the parton a is 
probed at high values of xi, it is legitimate to describe it with the collinear 
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parton density fa/Ai^i^l'''^)- Iii oui study we used one of the standard pdfs 
CTEQ6mE [3j. On the other hand, for the incoming gluon, which is probed 
at small X2, we must use the unintegrated gluon density (pgfB{x2, k'^), which 
in addition depends on that gluon's transverse momentum. The sum in 
Eq. ([2]) runs over the following sub-processes: qg — )• qg, gg — t- qq and 
99 99- 

3. Unintegrated gluon distribution in proton and nucleus 

In our study, we used the unintegrated gluon density from the unified 
BK/DGLAP framework [JHS]. This approach provides a well behaved, non- 
linear gluon distribution inside the proton and it incorporates the main 
sources of higher order effects. We generalized this framework to the case 
of a nucleus by assuming the Wood-Saxon nuclear density profile with the 
nucleus radius given by -Ra = RA^^^ where R is the proton radius and A 
is the mass number. In the limit A ^ 1 the result for the proton is recov- 
ered. The corresponding evolution equation for the distribution of gluons 
per nucleon in the nucleus A reads 



a.iVc f'dz r^dl^ f PHI,l'),Oi^ - I') - k'm,k^) , k^m,k') ] 
Jc, z P 1 \P-k^\ |4Z4 + A:4|3 J 




(3) 

where we used the simplified notation: (p{x,k'^) = (t>g/A{x,k'^) and ag = 
as(fe^). The second line in Eq. ([3|) corresponds to the BFKL equation with 
the kinematic constraint [1] (introduced by the theta function). The third 
line supplements DGLAP-type corrections and the last line corresponds 
to the nonlinear term whose strength is enhanced for the nucleus by the 
factor A^/'^. 

We started off by taking Eq. ([3]), with an appropriate initial condition, 
and fitted it to the most recent combined HERA data [7], which allowed us 
to fix the free parameters of the framework (proton radius R and three other 
coming from the parametrization of the initial condition). We performed 
the fit in the kinematical range of x < 0.01 and the full range of and 
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Fig. 2. Ratio of d-Au/p-p coincidence probabilities C(A(/)) for the forward di- 
hadron production at RHIC as a function of the azimuthal distance between the 
particles. The band corresponds to our prediction with the uncertainty related to 
unknown yield of uncorrelated dihadron production. For full details see [5] . 

obtained a very good description of data, with x^/ndof = 1.73. For detailed 
discussion we refer to [8]. 

4. Signatures of saturation in dijet production in p-A collisions 

Once the framework for the unintegrated gluon density had been fully 
specified, we used it to study saturation effects in the dijet production in 
p-A collisions. One observable which is very well suited for studying these 
effects is the distribution of dijet azimuthal distance Acp. In the region 
Acp ~ IT, which corresponds to the two jets being produced almost back-to- 
back, the gluon density is probed at low kt where the nonlinear effects are 
strong and the gluon distribution is expected to be significantly suppressed. 

Even though our main focus was on the central- forward dijet production 
at the LHC, the framework described above allowed us also to make an 
estimate of the suppression of dihadron production in d-Au collisions at 
RHIC. Fig. [2] shows the ratio of the normalized vr'' yields in the d-Au and 
the p-p collisions from STAR [9j. There, we also show our prediction which 
follows directly from using Eq. ([2]) together with the unintegrated gluon 
distribution from Eq. ([3]) with all the parameters fixed by HERA data and 
the mass number set to A=196 (for detailed procedure see [8j). As we see, 
our prediction correctly reproduces the suppression observed at RHIC in 
the region A(p ~ vr, which shows that our theoretical framework captures 
the essential physics of this class of processes. 

We then moved to the central-forward dijet production in the p-Pb col- 
lisions at the LHC. In Fig. Owe show the corresponding cross section as a 
function of the azimuthal distance between the jets for two energies of the 
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Fig. 3. Differential cross sections for central-forward dijet production at ^/s = 
5 TeV and 8.8 TeV as functions of azimuthal distance between the jets A(/) for the 
case of p-p and p-Pb collisions and two different cuts on jets' pt- 

p-Pb collisions, i.e. the current ^/s = 5 TeV and the nominal ^/s = 8.8 TeV. 
We also used two different jet pt cuts 15 and 25 GeV. 

We see that the A(/> distribution in the peak region is suppressed by the 
factor two for the case of the p-Pb collision with respect to p-p, and the 
effect extends to lower values of /S.(f) as we lower the jet pt threshold. This is 
a consequence of gluon saturation which is stronger in the Pb nucleus then in 
the proton and leads to suppression of gluon density at low kt, which is the 
region probed by the dijet configurations with ~ vr . As expected, the 
total yields increase with energy and decrease with pt cut but the relative 
difference between the p-p and p-Pb case remains similar. 

The distributions shown in Fig. [3] could be subjet to further corrections 
coming from the Sudakov and parton shower effects. They are however 
expected to act in a similar way for the proton and for the heavy ion collision 
since they affect the hard scattering. We would like to emphasize that the 
suppression observed in Fig. [3] comes from the initial state parton density 
saturation and therefore the relative difference between the cases of p-p and 
p-Pb will persist even if the very small region near = vr may profit from 
further refinements. 

5. Conclusions 

We presented predictions for the azimuthal dijet decorrelations in the 
p-p and p-A collisions in the framework of the high energy factorization 
with the unintegrated gluon density given by the nonlinear QCD evolution 
equation. This approach is unique as it allows one to study hard final states 
in the framework with saturation. 
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We validated our approach by using HERA F2 data and we checked 
that it correctly estimates the suppression of the away side peak observed 
by STAR. For the LHC, we found that the saturation in the Pb nucleus 
has a potential to manifest itself as a factor two suppression of the central- 
forward jet decorrelation in the region of the azimuthal distance between 
the jets A0 ~ vr. We argued that this relative difference should be largely 
insensitive to the final state effects. 
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